Critical slowing down exponents in quenched disordered spin models for structural 
glasses: Random Orthogonal and related models 



o 

(N 



F. Caltagirone^'^, U. Ferrari^'^, L. Leuzzi^'^, G. Parisi^'^'"^ and T. Rizzo"^'^ 

^ Dip. Fisica, Universitd "Sapienza", Piazzale A. Mora 2, 1-00185, Rome, Italy 

^ IPCF-CNR, UOS Rome, Universitd. "Sapienza", Piazzale A. Mom 2, 1-00185, Rome, Italy 

^ INFN, Piazzale A. Moro 2, 00185, Rome, Italy 

(Dated: January 30, 2013) 

An important prediction of Mode-Coupling- Theory (MCT) is the relationship between the power- 
law decay exponents in the /? regime. In the original structural glass context this relationship follows 
from the MCT equations that are obtained making rather uncontrolled approximations and A has 
to be treated like a tunable parameter. It is known that a certain class of mean-field spin-glass 
models is exactly described by MCT equations. In this context, the physical meaning of the so 
called parameter exponent A has recently been unveiled, giving a method to compute it exactly in 
a static framework. In this paper we exploit this new technique to compute the critical slowing 
down exponents in a class of mean-field Ising spin-glass models including, as special cases, the 
Sherrington-Kirkpatrick model, the p-spin model and the Random Orthogonal model. 
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I. INTRODUCTION AND FRAMEWORK 

It is well known that mean-field spin-glass models have 
a low temperature phase in which the Replica Symmetry 
is broken, with a breaking pattern that depends on the 
specific model. The models displaying a static discon- 
tinuous transition which are consistently described by a 
finite number of breakings are characterized by critical 
slowing down and a dynamical transition at a tempera- 
ture higher than the static one. 

They share some relevant properties of structural 
glasses [1-4], more specifically, the dynamical equations 
are exactly equivalent to those predicted by the Mode 
Coupling Theory (MCT) above the mode coupling tem- 
perature Tract where ergodicity breaking occurs. 
The time autocorrelation function in the high tempera- 
ture phase displays a fast decay to a plateau and then 
a second relaxation to equilibrium. Approaching the dy- 
namical transition temperature (called Td in the spin- 
glass context) the length of the plateau grows progres- 
sively until it diverges exactly at Td, where the system re- 
mains stuck forever in one of the most excited metastable 
states in a complex free energy landscape. According to 
MCT the approach to the plateau and the decay from it 
are both characterized by a power-law behaviour, respec- 
tively 



C{t) ~ gd + ct-" 



C{t) c^qa- c't" 



(1) 



(2) 



where qd is the height of the plateau and the two expo- 
nents satisfy the following relation that is exact in the 
framework of MCT (see for example [5]) 
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A 



(3) 



parameter A and, consequently, the exponents a and h 
have been computed exactly only for the spherical p-spin 
model [7] since the dynamical equations are particularly 
simple and correspond to the so called schematic MCT 
models. 

In most of the cases it is, instead, simply considered a 
tunable parameter, generically connected to the static 
structure function at Td through an often explicitly un- 
known functional [8]. In the case of continuous tran- 
sitions, instead, there is no dynamic arrest preceeding 
the static transition, the time correlation function does 
not display the two step relaxation and, consequently, 
no exponent h is defined. At the thermodynamic transi- 
tion, for long times, the correlation decays to the equi- 
librium value qEA with a power-law of the kind C{t) ~ 
qEA + ct^"" ■ The equilibrium order parameter qEA is zero 
at the transition in absence of magnetic field. 

It has been recently pointed out [9] that there exists 
a connection between the exponent parameter and the 
static Gibbs free energy, which allows to compute A in 
a completely thermodynamic framework, even in cases 
which go beyond schematic MCT. In the following we 
will briefiy summarize the method. 

Given a fully-connected model it is possible to com- 
pute the Gibbs free energy T{Q) as a function of the or- 
der parameter that, in the case of a spin-glass transition 
is the well known overlap matrix Q. The value of the or- 
der parameter can be determined through a saddle point 
calculation and T{Q) can then be expanded around this 
solution. For our "dynamic" purposes, the expansion has 
to be performed around a replica symmetric saddle point 
solution Qfff = q^^ . This gives raise to eight different 
kinds of third order terms, but only two of them will be 
relevant, namely: 



This relation has been proven to be robust under higher 
order corrections to standard MCT [6]. The exponent 
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In the case of discontinuous transitions it can be shown [9] 
that the following relation holds at the dynamical tran- 
sition, giving the connection between the dynamical ex- 
ponents a and h and the static coefficients, namely 



A = 



W2{Td) 

Wi{Td) 



(6) 



where A is given in Eq. (3) the expansion of the Gibbs 
free energy has to be performed around the value of the 
overlap yelding the height of the plateau at the dynamical 
transition, q^. 

Since the coefficients have to be computed at the dy- 
namical transition, where quantities at infinite time do 
not relax to their equilibrium (thermodynamic) value 
but remain stuck at their value inside the most excited 
metastable states, the averages should then be computed 
inside a single state. This corresponds to taking a 1-RSB 
Ansatz with breaking parameter tti — > 1 or, cquivalently, 
a RS Ansatz with the number of replicas n — > 1 [10-13]. 
In this paper we will use the second strategy which is 
technically much simpler than the first one, therefore we 
cannot treat the case of a dynamical transition in pres- 
ence of a magnetic field, since the mutual overlap (go) 
between states is non-zero and the 1-RSB (m -^ 1)/RS 
(n — > 1) equivalence docs not hold. 

On the other hand, for continuous transitions a relation 
between the exponent a and the two coefficients wi and 
■W2 analogous to Eq. (6) holds at the static point: 



A 



W2{Ts) 
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In this case, since the continuous static transition coin- 
cides with the dynamical one (e.g. in the SK model), 
the dynamical quantities at infinite time relax to their 
static value [14] and the averages can be computed in a 
replica symmetric Ansatz taking finally the limit n -^ 0. 
For this reason, if the transition is continuous, the RS 
Ansatz will be sufficient to treat the case in presence of 
a magnetic field. 

In order to compute the two coefficients wi and W2 one 
must determine the expression of the Gibbs free energy 
as a function of the overlap and then expand it to third 
order around the RS thermodynamic value q. In fully 
connected models, introducing a replicated external field 
e, the free energy reads 

f{£)^-J-^lnJdQe^pNiS[Q]+TTeQ) (8) 

which, for iV — > oo, can be evaluated at the saddle point 

/(£) = - ^extrg {S[Q] + Tr eQ) (9) 



We can immediately notice that the equation above ex- 
actly defines /(e) as the Anti Legendre Transform of the 
effective action 



f{e)= ALT S[Q] 



(10) 



and, again, by definition the Gibbs free energy T(Q) is 
the Legendre Transform of /(e), yelding 

r(0) = LT/(e) = LT (ALT5[Q]) = S[Q] (11) 

This implies that the functional form of the Gibbs free 
energy is equal to the one of the effective action. In 
fully connected models, we can then directly expand the 
latter. The general form of the third order term in the 
free energy is 



5(3) = 
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dQabdQcddQef 

Since a y^ b, c y^ d and e ^ f and the coefficients W 
are computed in RS Ansatz, we can have eight different 
vertices. It can be shown [15, 16] that, restricting the 
variations to the replicon subspace, i.e. the subspace 
where ^j, ^Qab = 0, one obtains the following expression 
containing only the two interesting coefficients wi and 
W2: 
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that follows quite straightforwardly from equation (12) 
applying the replicon constraint to the variations. 
In this paper we apply this technique to study the crit- 
ical slowing down of a general model of mean-field Ising 
sping-glass which includes, as particular cases, the SK 
model, the p-spin model and the Random Orthogonal 
model (ROM). The outline of the paper is the following: 
in section II wc introduce the general model, in section 
III we give the details of the computation of the param- 
eter exponent A for the general case and briefly present 
the result for the SK model and p-spin model. In sec- 
tion IV we compute A for the ROM model and compare 
our exact result with numerical simulations. Finally, in 
section V we give our conclusions and remarks. 



II. THE GENERAL MODEL 

In this section wc will consider a class of mean-flcld 
models with Hamiltonian 



n = -X!"^»J^»'^J"E 
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where Ui are N Ising spins. The 2-body interaction ma- 
trix is constructed in the fohowing way [17]: 



J^O^SO, 



(16) 



where O is a random 0{N) matrix chosen with the Haar 
measure. On the other hand, S is a diagonal matrix with 
elements independently chosen from a distribution p{^). 
In order to ensure the existence of the thermodynamic 
limit, the support of p(^) must be finite and independent 
of N. The p-body interactions K^ are i.i.d. gaussian 
variables with variance 
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for some real valued function K. As shown in [17-19] for 
this class of mean field spin-glasses, the general form of 
the replicated free energy is: 



-n/3/ = extrQ,A5[Q,A] 



(19) 



with 



5[0, A] - ^Tr G(/30) + y ^ R{Qa^ 

ah 
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where G : Mnxn ^ Mnxn is a (in general rather compli- 
cated) function in the space of n x n matrices, formally 
defined through its power series around zero. The par- 
ticular form of G depends on the choice of the eigenvalue 
distribution p{^). In the following we will consider mainly 
two cases: Wigner law and bimodal. 
Given this effective action, the saddle point equations in 
A and Q respectively read 



Qab = {{<^aO'b)) 
Kb = l3[G'{PQ)]ab 



P^R'iQab) 



(21) 



where the average ((•)) is computed with the measure 



>V(A,f7) = 






(22) 



In Replica Symmetric Ansatz {Qab =^ 1 , Kb = -^ for 
a 7^ 6 and Qaa = Qd , Ka = Ad), Eq.s (21) become 



(23) 



P^r.n 



A = ^G'{P{l + {n-l)q))-G'mi-q))]+p^R'{q) 
n 
where m = tanh(z) and the average (•) is computed with 

the measure 



/i(A) =: e 2A cosh"(z) 
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In the next section we study in detail the (dynamical) 
critical behaviour of this class of models and we show 
how to compute the critical slowing down exponents. 



III. COMPUTATION OF THE MCT 
EXPONENTS 



As explained in the Introduction, in order to compute 
the parameter exponent A we have to expand the effective 
action to third order in Q and then restrict the variations 
to the replicon subspace, obtaining straightforwardly the 
two coefficients Wi and W2- 

In the present case, the effective action contains the 
auxiliary field A which will be eliminated making use of 
the saddle point equation (21). The expansion of Eq. 
(20) to third order gives: 



5S[Q,K] ~;^Tr[/3G"(/3Q'5^),5g + J/32G(2)(/30^^)<5Qr5Q+ ^/33g(3)(/3Q^^)(5Q(5Q5Q] 
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A comment is needed for the first line of cq. (25): the 
"scalar like" Taylor expansion of a matrix functional 
/(M), / : Mnxn — > Mnxn around some Mq (different 
from the null matrix), is correct only if [Mq,SM] = 0. 

In the present case Q^^ is replica symmetric while 
6Q is, in principle, simply symmetric. The commutation 
condition for a RS matrix with a symmetric matrix reads: 



'^^^ 5Qcb = y^ SQac ya,b 









(26) 



that is satisfied in any subspace orthogonal to the anoma- 
lous one (see [15]), i.e. both in the longitudinal and in 
the rcplicon sector [27]. Equating to zero the first order 
of Eq. (25) we obtain the saddle point equations (23). 

Considering that the variations SQ and SA are in the 
rcplicon subspace the second order term simplifies as fol- 
lows 



[/?^(.9f-ff(^')+/3V^']E'5Qa. 
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where, here and in the following formulas we define the 
four constants (two diagonal and two off-diagonal): 
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For the system to be critical, the rcplicon eigenvalue must 
vanish and, consequently, the Hessian determinant must 
be zero. 

Imposing this condition we get the following equality 



{{l-m'Y) 
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[P'^igT -a^'^'^) + P'^T^'^^] 



(30) 



Eq. (30) together with Eq.s (23) gives the criticality con- 
dition, which locates the dynamical or static transition 
point depending on the value (n = 1,0) of the replica 
number. 

Now we want to eliminate the auxiliary field. The A - 
saddle point equation (21) reads (up to second order) 



cd 

+ I E ((TWV'^aV-f))c'JA,<idAe/ 



(31) 



Exploiting the property of the replicon subspace and the 
criticality condition we can write the variation in the fol- 
lowing way: 

5Qab = {{l-m'f)5Kab 

+ {{l-m'f)Y,5Kac5Kb (32) 

C 

Inverting the equation we obtain 

<5A„, = [/32(gf -.g(2)) + /3V^)]jg,, 

-[,3^(5f-.g(^)) + /3V2']-^x 
\{l-m'f)Y,5Qac5Qcb+ 



(33) 



Now we define 



2{m\l-m'f)5Ql,) 
C2 = 2(m2(l-m2)2) 



(34) 



and plug the constraint (33) into (25), obtaining three 
different contributions to the third order in JQ, namely 
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Summing all the third order contributions in Eq. (25), 
we eventually find 



3!' 
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Substituting Eq. (28) one immediately obtains a general 
expression for the exponent parameter, that is the main 
result of this paper: 
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(40) 



R{x) = ia;2 so that R^'^\x) = 1 and R'^^^x) = 0. Then 
we get immediately the expression for the exponent pa- 
rameter: 



_C2 _ 2(m2(l-m2)2) 
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Our result exactly coincides with the one obtained by 
Sompolinsky and Zippclius [14] in a purely dynamical 
framework. 



SK model on the dAT line 



B. Multi p-spin Ising model 



The Sherrington - Kirkpatrick model [20] is described 
by the Hamiltonian 



n 



- ^ X! ■^^3 ^^'^^ ^ ^^'^^ 



(41) 



where the couplings are i.i.d. random variables dis- 
tributed according to a gaussian with zero mean and 
variance 1/A^. This model belongs to the class defined 
above, with R — x^ /2 and G = or, conversely, R — 
and G = x'^/2 [20], except for the presence of the mag- 
netic field term. We will see in a while that this affects 
the result in a very simple way. 

It is well known that in the SK model there exists a line of 
instability of the replica symmetric solution in the j3 — H 
plane, the de Almeida - Thouless (dAT) line [21], where 
the so called replicon eigenvalue of the stability matrix 
vanishes. In this section we want to compute the decay 
exponent of the time correlation function along this line. 
In order to get the result we first have to find solutions 
which satisfy simultaneously the saddle point and the 
dAT equations, respectively 



dii{z) tanh^(/3ygz + (iH) 
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Notice that the expression (39) is completely general 
and holds for every model belonging to this class, while 
the details of the model enter in the specific form of the 
functions G and R. In the case of continuous transi- 
tions, the presence of the magnetic field, only modifies 
the definition of the parameter m in equation (39) which 
becomes: 



Starting from an Hamiltonian of the kind of Eq. (15) 
without the first term, leads to a generalized version of 
the p-spin model [22], in wich many multibody interac- 
tion terms arc considered, depending on the actual form 
of the function R. As shown in refs [23-25] , in these mod- 
els, the thermodynamic properties, the critical dynamics 
and replica symmetry breaking structure depend on the 
relative strength of the coupling terms (the coefficients 
of the expansion of R). Indeed, in order to treat a par- 
ticular case, before applying our technique, one should 
understand the behavior of the corresponding model. 

The simple p-spin Ising model is characterized by 



Rix) 



where a o 7^ 1 affects only the variance of 



the couplings distribution and indeed rescales the tem- 
perature. For p > 2, in absence of any external mag- 
netic field, the model displays a standard RFOT transi- 
tion from the paramagnetic RS phase to the 1-RSB spin- 
glass and, at a lower temperature, a second transition to 
a FRSB spin-glass [22]. Focusing on the first transition, 
we can compute the critical dynamic exponents in the 
present general framework while a specific analysis was 
presented by some of us in [16]. In order to recover the 
exact same model, we have to set R{x) ~ xP/2, which 
reduces Eq.(23) and (30) to: 



A = ^,-^ 
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and Eq. (39) to: 
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as it was found in [16]. 



IV. RANDOM ORTHOGONAL MODEL 



m = tanh(z + f3H) 



(44) 



without changing the formal expression for the coeffi- 
cients. Therefore, plugging into the effective action (20) 



The Random Orthogonal Model (ROM) [17, 18] is ob- 
tained with the choice R = and 



pi0^c,S{^-l) + il-a)5(i + l) 



(48) 





FIG. 1: The static and dynamical critical temperature as a 

function of the parameter a. 

They both diverge for a — > and are zero at a = 1. 



FIG. 2: The static and dynamical critical overlap as a function 
of the parameter a. 



It displays a random first order transition regardless of 
the value of the tunable parameter a. The case with 
a = 1/2 has been extensively sudied in [18] while the 
general case was treated in [17]. 

They show that, as a consequence of the choice of the 
eigenvalue distribution (48), the function G appearing in 
the effective action reads [28]: 

2G{x) = [1 + 4.T(2a - 1 + x)]^^^ 

+ (2a - 1) ln{[l + 4:x{2a - 1 + .t)]^/^ + 2a; + 2a - 1} 

- ln{[l + 4a;(2a - 1 + x)]^/^ + 1 + 2(2a - l)x} 

- (2a - 1) ln(2a) - 1 + ln(2) 

(49) 

Using Eq.s (21) we can determine the transition temper- 
ature Td and the dynamical overlap qd which are shown 
in Figs. 1 and 2 and coincide with those found in Ref. 

[17]. 

Once the critical point is obtained as a function of 
a, using formula (39) specialized to the ROM case, we 
obtain the value of the exponent parameter A (a) and of 
the critical slowing down exponents a(a) and 6(a) which 
are shown in Fig. 3. We find numerically that for a — > 1 
the exponent parameter goes to A = | as in the Ising p- 
spin model for p — > oo [16]. On the other hand for a — > 
we find A = 2 ^s ii^ the Ising p-spin model for p — > 2. 

In particular, for a = 13/32 we have b = 0.628. We 
now use this value to compare with numerical simula- 
tions. 



A. Comparison with Monte Carlo data 

There are recent numerical simulations by Sarlat et al. 
[26] on the fully connected ROM which give an estimate 




FIG. 3: Solid line: the exponent parameter A. Dashed line: 
the exponent a. Dot-dashed line: the exponent b. 



for the MCT exponents a, b and 7. They choose a — 

1| ~ 0.4 in order to have higher transition temperatures 

and a good separation between the static and dynamical 

critical temperature. 

Their direct estimate of the exponent b is 0.62, while their 

direct estimate of 7 is 2.1 which, through the exact MCT 

relations: 



T^{i + b) _ r^(i 
f(i 
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yields bucT — 0.75. 

Our exact computation yields instead btt — 0.628, which 
suggests that the best estimate of the exponent b in [26] 
is the direct one, that is very close to the actual value. 



V. SUMMARY AND CONCLUSIONS 

In the present work we have introdueed a general fuUy- 
connected model for Ising spins, which combines an or- 
thogonal two body interaction with a set of p-body in- 
teractions. 

Exploiting a technique that has been recently intro- 
duced [9], based on the equivalence between statics and 
long time dynamics, we have been able to find an an- 
alytic expression for the exponent parameter A, in a 
purely static framework. As particular cases of the gen- 
eral model we have studied the Sherrington-Kirkpatrick 
model along the de Almeida- Thouless line, the p-spin 
model and the Random Orthogonal model. 

For the SK model we find the same result found by 
Sompolinsky and Zippelius in Ref. [14]. 

For the p-spin we recover, as a byproduct of the gen- 
eral model, the results given in detail in ref. [16]. 



We have studied the critical behaviour of the param- 
eteric class of Random Orthogonal models at arbitrary 
values of the constant a G [0,1], that determines the 
distribution of the eigenvalues of the interaction matrix. 
The exponent parameter and the two MCT exponents 
have been determined analitically for any a and in par- 
ticular we have looked at a — 13/32 in order to make 
a comparison with existing numerical simulations [26]. 
Our exact result is in very good agreement with the one 
obtained in the Monte Carlo study, through a direct es- 
timate of b (late /3 regime). 

On the other hand, a direct estimate of the exponent 7 
gives a result that is quite far from what we found here, 
suggesting that the strong finite-size corrections affect 
the value of 7 much more than b. Numerical interpo- 
lations at criticality are very sensitive for glassy models 
and the corresponding estimates can strongly suffer of 
this drawback. 



[9 

[10 

[11 
[12 

[13 
[14 
[15 
[16 



T. Kirkpatrick and P. Wolynes, Phys. Rev. B 36, 8552 [17 
(1987). 

T. Kirkpatrick and D. Thirumalai, Phys. Rev. B 36, 5388 [18 
(1987). 

T. Kirkpatrick and D. Thirumalai, Phys. Rev. Lett. 58, [19 

2091 (1987). [20 

T. R. Kirkpatrick and P. G. Wolynes, Phys. Rev. A 35, 
3072 (1987). [21 

W. Gotze, Complex Dynamics of Glass-Forming Liq- 
uids: A Mode-Coupling Theory (Oxford University Press, [22 
2009). [23 

A. Andrcanov, B. G., and B. J.-P., EPL 88, 16001 (2009). 
A. Crisanti, H. Horner, and H. Sommers, Z. Phys. B 92, [24 

257 (1993). 

F. Weysser, A. M. Puertas, M. Fuchs, and T. Voigtmann, [25 
Phys. Rev. E 82, 011504 (2010). 

F. Caltagirone, U. Ferrari, L. Leuzzi, G. Parisi, F. Ricci- [26 
Tersenghi, and T. Rizzo, Phys. Rev. Lett. 108, 085702 
(2012). [27 

A. Crisanti, Nucl. Phys. B 796, 425 (2008). 
R. Monasson, Phys. Rev. Lett. 75, 2847 (1995). 
S. Franz and G. Parisi, J. Phys. I (France) 5, 1401 (1995). 
S. Franz, G. Parisi, F. Ricci- Tersenghi, and T. Rizzo, 
European Physical Journal E 34, 1 (2011). 

H. Sompolinsky and A. Zippelius, Phys. Rev. B 25, 6860 [28] 
(1982). 

T. Temesvari, C. De Dominicis, and I. Pimentel, The 
European Physical Journal B 25, 361 (2002). 
U. Ferrari, L. Leuzzi, G. Parisi, and T. Rizzo, 
arXiv;1202.4168vl (2012). 



R. Cherrier, D. S. Dean, and A. Lefevre, Phys. Rev. E 

67, 046112 (2003). 

E. Marinari, G. Parisi, and F. Ritort, Journal of Physics 

A: Mathematical and General 27, 7647 (1994). 

A. Crisanti and H. Sommers, Z. Phys. B 87, 341 (1992). 

D. Sherrington and T. Kirkpatrick, Phys. Rev. Lett 35, 
1792 (1975). 

J. R. L. de Almeida and D. J. Thouless, J. Phys. A 11, 
983 (1978). 

E. Gardner, Nucl. Phys. B 257, 747 (1985). 

A. Crisanti and L. Leuzzi, Phys. Rev. Lett. 93, 217203 
(2004). 

A. Crisanti and L. Leuzzi, Phys. Rev. B 76, 184417 
(2007). 

A. Crisanti, L. Leuzzi, and M. Paoluzzi, Eur. Phys. J. E 
34, 98 (2011). 

T. Sarlat, A. Billoire, G. Biroli, and J. P. Bouchaud, J. 
Stat. Mech. (2009). 

Any element of the anomalous subspace can be described 
by a one index field, i.e. by a vector tp" restricted to the 
condition ^^ ■tl)°'—0. A generic anomalous field can then 
be written as i>'f = 1/2(V'" + Tp'^). 

Any element of the longitudinal subspace can be de- 
scribed by a scalar 4'°'^ = Tp. 

Note that the corresponding formula in [17] contains a 
typing mistake: ln(2) is taken with the negative sign, 
which would give a negative entropy at infinite tempera- 
ture. 



